Mon. Not. R. Astron. Soc. 000, [THS] (2008) 



Printed 19 November 2008 



(MN I*TeX style file v2.2) 



Impact of sub-solar metallicities on photometric redshifts 



Ralf Kotulla* and Uta Fritze 

Centre for Astrophysics Research, University of Hertfordshire, 



Accepted XXX. Received XXX; in original form 2008 May 29 



College Lane, Hatfield ALIO 9AB, United Kingdom 



ABSTRACT 

With the advent of deep photometric surveys the use of photometric redshifts. obtained 
with a variety of techniques, has become more and more widespread. Giving access to 
galaxies with a wide range of luminosities out to high redshifts, these surveys include 
many faint galaxies with significantly sub-solar metallicities. 

We use our chemically consistent galaxy evolutionary synthesis code GALEV to 
produce a large grid of template Spectral Energy Distributions (SEDs) for galaxies of 
spectral types E and Sa through Sd - one accounting in a chemically consistent way 
for the increasing initial metallicities of successive stellar generations, the other one 
for exclusively solar metallicities - for comparison. 

We use our new photometric redshift code GAZELLE based on the comparison of 
observed and model SEDs. Comparing the photometric redshifts obtained using solar- 
metallicity templates when working on a catalogue of artificially created chemically 
consistent SEDs. typical for low-metallicity local late- type galaxies and for intrinsically 
low-luminosity, and hence low-metallicity, galaxies in the high-redshift universe, we 
find a significant bias resulting from this metallicity mismatch. This bias consists in 
a systematic underestimate of the photometric redshift by typically Az « 0.1 .. .0.2 
until z K, 1.2, depending on galaxy type, of distant, faint and low-metallicity galaxies 
if analysed with solar- metallicity templates. 

Key words: galaxies: evolution - galaxies: abundances - galaxies: high-redshift - 
galaxies: distances and redshifts 



, 1 INTRODUCTION 

. Most of today's studies of large samples of high-redshift 
' galaxies rely on photometric redshifts comparing observed 
spectral energy distributions (SEDs), i.e. magnitudes and 
colours in multiple filters, to a set of templates. For that 
reason the choice of the right set of template is crucial to 
determine accurate and unbiased photometric redshifts. The 
most widely use d templates are either observed local tem- 
plates, e.g. from lColeman et al.l (|l980l ), that do not includ- 
ing any evolutionary correction, or templates generated with 
evolutionary synthesis models which in most cases use fixed 
solar metallicities. 

During the last years, evidence has accumulated that 
galaxies are not made out of stars of one metallicity, but 
show a wide range from very metal-poor to more metal- 
rich stars. This holds true not only for ou r Milky Way 
(|Rocha-Pinto fc Maciell Il998l : lAk et al.ll2007l). but also for 
external galaxies like e.g. the LMC ( Cole et al]|2000i) and gi- 
ant ellipticals like NGC5128 {— Centaurus A, Harris et al.l 
1 19991 : [Harris fc Harrisllioool ). 

Studying samples of local star forming galaxies 
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ISkillman et al.l (|l989l ) showed a trend of decreasing average 
metallicity of a galaxy with decreasing luminosity, used as 
an indicator of its mass, and spanning more than 12 magni- 
tudes i n luminosity. Larger samples, compiled e . g. fro m the 
SDSS (|Tremonti et all |2004| : iKewlev fc EllisonI |2008| ) con- 
firmed this m ass-met allicity relation and extended it to even 
lower masses (|Lee et al. 2006). 

This is particularly important with respect to high 
redshift galaxies in the early universe, since those galax- 
ies did not yet have the time to produce enough stars to 
enrich their IS M to high metallic ities. Studies of Lyman 
Break G alaxies (iPettini et al.ll2p01^ . Damped Lyman-a Ab- 
sorbers dProchaska et alT 20031 ) and Gamma Ray Bursts 
(|Prochaska et al.ll2007F all show that galaxies get progres- 
sively more metal-poor if they are observed at high red- 
shift. Furthermore dwarf galaxies, such as the LMC, that 
are metal-poor in the local universe, are observable out to 
considerable red shifts in deep imaging surveys. Note that 
lErb et al] (|2006l ) found a mass-metallicity relation for galax- 
ies at redshifts of 2 ~ 2, confirming that the local trend was 
already established in the early universe. 

For that reason evolutionary synthesis models that 
take the chemical enrichment of successive stellar genera- 
tions into account are principally superior to more simpli- 
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fied models with fixed metallicity. In iBicker et al.l (|2004l ) 
we showed that with appropriate star formation histories 
(SFHs) our chemically consistent GALEV models agree well 
with a wealth of observed properties for local and high- 
redshift galaxies and showed how the presently observed 
ste llar metallicity distributions in galaxies have evolved. 
In iBicker fc Fritzd (|2005l ) we demonstrated the impact on 
non-solar metallicities on the determination of star for- 
mation rates ( SFRs ) from emission lines and UV-fluxes. 
iKodama et~al] (|l999l ) have shown that metallicities can have 
a significant impact on observed colours of high-redshift 
galaxies; the impact of metallicity on restframe colours was 
studied by other authors. However, the impact of those gen- 
erally bluer colours on photometric redshifts has not been 
studied so far. In the present Letter we quantify the impact 
of neglecting those sub-solar metallicities on photometric 
redshifts. 



2 CREATION OF TEMPLATE SEDS 
2.1 Input models 

To study the chemical enrichment history of the common 
spectral galaxy types E and Sa through Sd we used our 
chemically consistent galaxy evolution code GALEV. As- 
suming a closed-box model, GALEV allows us to compute 
the chemical enrichment of a galaxy's gas-reservoir from the 
yields of dying stars. We use isochrones from the Padova- 
group (|BertelIi et al.|[r994l ) with metallicities ranging from 
[Fe/H] = -1 7 to [Fe/H] = +0.3 and a Salpeter-IMF 
iSalpeteij [l955l ') with mass-limits of 0.10 Mq and 100 Mq. 
Note that a different choice of the IMF, e.g. Kroupa or 
Chabrier, does not affect the results obtained below. The 
spectral galaxy types are characterized by an exponentially 
declining SFR for the E-model, SFRs proportional to the 
available gas-mass for the Sa-Sc models (with factors of pro- 
portionality decreasing towards later types), and a constant 
SFR scenario for the Sd. These SFHs were shown to pro- 
vide a good mat ch to today's ga l axy t emplates, e.g. from 
iKennicutd (| 19921 ') in IBicker etall (|2004 ). To derive the ef- 
fects on photometric redshifts we also ran all those mod- 
els again, this time fixing the metallicity to the solar value 
[Fe/H] = 0.0. Our GALEV models also include line and 
continuum emission from ionized gas, with line strength ap- 
propriate for the particular gaseous metallicity at each time. 
The importance of emission lines for accurate photome tric 
redshifts has recently been shown bv lllbert et all (|2008l '). 

To derive the SEDs for comparison with observations, 
we redshifted the spectra at all timesteps to the redshift 
corresponding to the age of each galaxy spectrum, adopt- 
ing a concordance cosmology with Ho = 70kms^^ Mpc^'^, 
Qui ~ 0.30 and Q\ = 0.70. We assume that all galax- 
ies started forming stars at 2 = 8. Variations of this for- 
mation redshift wi th spectral type, as e.g. suggested by 
iNoeske et all (|2007l l?h for the latest spectral types has little 
effect on the present study since those show little evolution 
anyway. This means we consistently account for both evo- 
lutionary and cosmological corrections. Since some galaxies 
at high-reds hifts are found to contain significant amounts 
of dust fe.g. lSteidel et al.lll999l ). we convolve each spectrum 
with the dust attenuation curve of ICalzetti et al.f (|2000l ). 



We chose a range of extinctions from E(B — V) = 0.0 mag 
to E(B - V) = 0.5 mag in steps of AE(B - V) = 0.05 mag. 
We include th e effects of intergalactic absorption following 
iMadaul j 19951 ) and then convolve each spectrum with a set 
of typical filters, here taken to be the SDSS ugriz and the 
standard 2MASS JHK filters. Galaxies were normalized as 
to show by z=0 the average observed Johnso n B-band mag- 
nitudes for their respective types as given bv lSandage et all 
(|l985al lbl). We then added the bolometric distance modulus 
for each redshift. This results in a total of ~ 3200 template 
SEDs with smaller redshift intervals at lower z and wider 
sampling at high redshifts for each of our ten models (5 
types E, Sa-Sd, all chemically consistent and with metallic- 
ity fixed to solar for comparison). 

2.2 Addition of noise 

To simulate real observations, we added Gaussian noise ob- 
servational errors to each magnitude. The amount of scatter 
added, Ami, depends on the magnitude mi of the i-th filter: 

Ami = a -f b X exp(c x mi — d), 

with a = 0.03 mag describing calibration or zeropoint uncer- 
tainties, b — 3.75 and c = 0.75 being constants defining the 
shape of the curve, d depends on the depth of the under- 
lying observations, here chosen to correspond to 5a limit- 
ing AB magnitudes of (26,27,27,27,26,25,25,24) mag for the 
(u,g,r,i,z,J,H,K) filters. This procedure was repeated 100 x 
for each input SED, resulting in an artificially created in- 
put catalog of ~ 7 X 10^ galaxies for each model. For the 
following analysis we then derived median values and 1 a- 
uncertainties in bins of Az = 0.05. 

2.3 Determining photometric redshifts 

To derive the photometric redshifts we use our photomet- 
ric redshift code G AZELLE described in more details in 
a companion paper (|Kotulla fc Fritzell200j '). In principle it 
uses a x^-algorithm to compare fluxes derived from the ob- 
served SEDs with a range of template SEDs. The resulting 
X^-values are then transformed into normalized probabili- 
ties. Masses are derived by scaling the model SEDs as a 
whole to match the observed SED on average. To deter- 
mine Icr-uncertainties for redshifts and all dependent pa- 
rameters (masses, SFRs, metallicities, etc.) we derive the 
minimum and maximum values encountered while summing 
up normalized probabilities (from highest to lowest) until 
68 % have been reached. We restrict our template set to 
only undisturbed galaxies E, and Sa through Sd, since those 
are well calibrated against observed galaxy templates, and 
match observations in colours, spectra and metallicities (see 
IBicker et al. I I2OO4I : iKotuUa et al. I I2OO8I . for a detailed com- 
parison) . 

In the following, we will focus on the redshift determi- 
nation and the best-match x^-value. 



3 RESULTS 

3.1 Evolution of metallicity with redshift 

In Fig. [1] we present the metallicity evolution of the dif- 
ferent spectral galaxy types E, Sa, and Sd with decreasing 
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Figure 1. Chemical enriciimcnt histories for galaxies of differ- 
ent spectral types E (top panel), Sa (lower panel, upper curves) 
and Sd (lower panel, lower curves). The blue solid lines marks 
the gas-phase or ISM-metallicity, while the dashed lines represent 
luminosity-weighted stellar metallicities in different bandpass es. 
Black points mark observed metallicities l lZaritskv et al ][i99l of 
local galaxies. 



redshift. We show two different metallicity measures: The 
gas-phase or ISM metallicity and the luminosity-weighted 
stellar metallicity in a set of different rest-frame filters. 
The ISM metallicity is traditionally measured from emis- 
sion lines, while stellar metallicities are deriv ed from stel- 
lar absorption lines, as e. g. Lick indices (e.g. iTrager et al.l 
ll998l : ISchiav"on" et al I2OO6I ). requiring spectra of much higher 
signal-to-noise ratio. Our models yield metallicities at z = 
of Ze = Zq, Zsa = 1.5 Zq and Zsd = 0.25 Zq, 
in good agreement with observed metallicities, e.g. from 
IZaritskv et ai] (|l994h . 

The most important point with respect to this paper, 
however, is that only the E and Sa models reach enrichment 
levels comparable to solar metallicity. Later spiral types, i.e. 
the Sb-Sd models only reach significantly sub-solar metallic- 
ities after a Hubble time, rendering the assumption of solar 
metallicity independent of redshift and galaxy type invalid. 
While solar metallicity is a moderately good approximation 
(Z(t) > 0.5 Z0) for early-type galaxies (E to Sa) back to 
fairly young ages or high redshifts, it becomes less and less 
valid for later galaxy types, in particular at earlier times or 
higher redshifts. 

Low-metallicity stellar populations are brighter in the 
optical and UV (cf. Fig. [2]), have bluer colours (since their 
stars are hotter) and produce more ionizing photons com- 
pared to their equal-mass solar metallicity counterparts. 
This leads to higher emission line fiuxes an d hence an over- 
estim ation of SFRs by up to factors ^ 2 (|Bicker &: Fritzd 
I2OO5I ) if solar-metallicity calibrations are used. At the same 
time, their higher overall luminosities lead to overestima- 
tions of galaxy masses by factors up to ^ 5 and their bluer 



1 








CC " 

^solar 


0.1 










0.01 










0.001 












1000 


2000 


4000 


10000 20000 



wavelength (A) 

Figure 2. Spectrum of a 4Gyr old constant star formation rate 
model (Sd) calculated in the chemically consistent way (upper 
blue curve) and the fixed solar metallicity only way (lower red 
curve), both having identical masses. 

colors lead to an underestimation of their stellar population 
ages by factors up to ^ 2, unless their sub-solar metallicities 
are properly taken into account. 

3.2 Impact on photometric redshifts 

We ran three different sets of photometric redshift deter- 
minations, comparing (a) the solar metallicity SED catalog 
to solar metallicity SED templates, (b) the chemically con- 
sistent (CC) SED catalog to CC SED templates, and (c) 
analysing the CC SED catalog using solar metallicity SED 
templates. The outcomes of runs a) and b) are shown by the 
green and blue lines and symbols in Figures [3] and |31 every 
datapoint represent the median-value in bins of Az = 0.05 
in redshift. As expected, we find very small x^-values for the 
best match and excellent correspondence between true and 
photometric redshifts. The third run analyzing the CC cata- 
log with solar metallicity templates mimics the wide-spread 
analysis method for observation of low-metallicity galaxies in 
the early universe using close-to-solar metallicity SED tem- 
plates. Those can either be locally observed galaxies that 
naturally have higher metallicities than their high-redshift 
counterparts, training sets of galaxies with available spec- 
troscopic redshifts (and hence the brightest and with the 
mass-metallicity relation also most metal-rich galaxies at 
each redshift) or solar metallicity model templates. The re- 
sult are shown as red symbols in both figures. As expected, 
the best- match x^/DOF- values (where degrees-of- freedom 
(DOF) means the number of filters) for run (c) are sig- 
nificantly larger at almost all redshifts. The trend towards 
smaller values at higher redshifts can be understood as a 
consequence of photometric uncertainties increasing with de- 
creasing brightness and finally a decreasing number of filters 
due to dropouts and magnitudes falling below the detection 
limit. This in turn allows more fiexible matching by varying 
both shape, determined by galaxy type, redshift and extinc- 
tion, and normalisation, i.e. galaxy mass, of the template 
SED. 

Fig. Ushows the offsets between true and retrieved pho- 
tometric redshifts that result from the choice of templates 
not matching the observed metallicities. We also show the 
1(7 regions for each galaxy type as filled regions. It is obvious 
that even for the near-solar metallicity galaxy types E and 
Sa there are still significant offsets of Az = Zphot — Zspec ~ 
-0.1 (equivalent to a^, = Az/(l-|-z) > -0.05) until z 1.2. 
At higher redshifts z = 1.5. ..2. 8 we also find a bias but 
this is less prominent than at lower redshifts, in particu- 
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Figure 3. x^"Valuc per degree of freedom (DOF) of best- 
matching galaxy-template combination as function of redshift for 
three different galaxy types E (solid, red line), Sa (blue dash- 
dotted line) and Sd (green dashed line). In all three cases we 
used solar-metallicity templates and chemically consistent input 
galaxies. Dark and light red shaded regions mark the Icr and 3cr 
ranges for the E-type model. The blue shaded region marks the 
outcome of the matched template runs (solar vs. solar and CC 
vs. CC). 



lar compared to the increased scatter at those redshifts. At 
even higher redshifts z > 3 dropouts start to dominate the 
redshift determination. The reason for these biases is that al- 
though the metallicity is near solar for the early types E and 
Sa, the galaxy nevertheless contains a large fraction of lower 
metallicity stars (e.g. « 2/3 of the U-band flux of nearby 
elliptical galaxies is emitted by stars with [Fe/H] ^ —0.7; 
iBicker et al.ll20o3 ). As a general trend, the retrieved pho- 
tometric redshifts show a bias towards lower redshifts. This 
trend can be understood from the bluer SEDs of the CC- 
models, generated by the lower metallicity stars, that the 
photometric redshift code tries to compensate for by at- 
tributing lower redshifts to its Zq templates. This offset is, 
in particular at low redshifts z ^ 1, larger than the typical 
scatter of CTz ^ 0.03 found for large samples of photometric 
redshifts (e.g. iMobasher et al.|[2007h . 

Furthermore, the amount of dust reddening also plays a 
role in the following sense: essentially dust-free models with 
E(B— V) < 0.1 mag have larger best-match values and are 
more strongly biased towards too low redshifts than models 
with more dust, reaching a maximum for the Sd-type galaxy 
at 2 « 0.4 where Az = —0.2 or equivalently — 0.14 (see 
dahed lines in Fig. |4ll . This point is particularly crucial at 
z > 1 where we essentially observe the rest-frame UV. We 
here are biased towards UV-bright objects, i.e. those that 
are not hidden behind large amounts of dust. Photometric 
redshifts obtained by fitting solar-metallicity templates to 
those galaxies are therefore even more strongly biased to- 
wards too low redshifts than the median of all extinctions 
presented above. 

Observation al evidence for the bias described here can 
be found e.g. in lllbert et al ] (|2006l . Fig.3). There observed 
templates were used to derive photometric redshifts from a 
filter set similar to the one used here, and a underestimation 
until z 0.6 and in particular at Zphot = 0.3, Zgpcc ~ 0.4 was 
found. 
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Figure 4. Redshift offset Az = Zpjjot — ztrue for the elliptical 
(top panel) and spirals Sa, Sb, Sc, and Sd (lower panel). In the 
top panel, blue and green lines are for matching combinations, red 
symbols show the errors resulting from the use of solar metallicity 
templates for the analysis of lower metallicity galaxies. Each point 
represents the median value in a bin of width Az = 0.05. The 
dashed lines in the lower panels show the bias for E(B — V) ^ 
0.1 mag. 



4 CONCLUSIONS AND SUMMARY 

We used our chemically consistent galaxy evolutionary syn- 
thesis models GALEV to study the chemical enrichment his- 
tories of galaxies over a range of spectral types E and Sa 
through Sd. The E-type galaxy reaches enrichment levels of 
Z > 0.5 Z@ already at high redshifts z « 4 and remains 
almost from there on. Sa-type galaxies are significantly sub- 
solar ai z > 1.5, while later types such as Sd even after a 
Ifubble time only reach levels of 1/4 Zq. 

This fact, in combination with observational evidence 
for a wide range in stellar metallicities of our and nearby 
galaxies and the decreasing stellar metallicities in galaxies at 
higher redshifts casts doubt on widespread methods of using 
only solar-metallicity templates to derive photometric red- 
shifts from observed spectral energy distributions. We study 
the impact of the increasing importance of sub-solar metal- 
licity populations in high-redshift galaxies on photometric 
redshift determinations using our photometric redshift code 
GAZELLE on several large synthetic galaxy catalogs and 
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find a significant bias of Az ~ 0.1 for galaxies at z ^ 1.2 to- 
wards systematically underestimated photometric redshifts 
as a consequence of their bluer SEDs. 
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